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Abstract 

A calculation of the heavy quark potential at finite temperature at strong coupling based on the 
AdS/CFT correspondence is presented. The calculation relies on the method of complex string 
trajectories previously introduced in JH], and on the introduction of a modified renormalization 
subtraction. The obtained potential is smooth, negative definite for all quark-antiquark separa- 
tions, and develops an imaginary part for r > r c - 0.870/nT . At large separations the real part 
of the potential does not exhibit the exponential Debye falloff expected from perturbation theory 
and instead falls off as a power law, proportional to 1/r 4 . 



The heavy quark potential at finite temperature is a very important quantity in the study of 
Quark Gluon Plasma formation in Heavy Ion Collisions carried out at RHIC. Until recently it has 
been calculated either analytically at small coupling using perturbation theory, or numerically us- 
ing lattice simulations, hence in Euclidean time. However, the matter formed in RHIC collisions 
is rapidly evolving in time and, according to latest analyses, it is strongly coupled. Thus, cal- 
culational tools for strongly coupled, real-time QCD dynamics are needed. Such tools are now 
available via the Anti-de-Sitter space/conformal field theory (AdS/CFT) correspondence, albeit 
for a theory different to QCD, namely N = 4 supersymmetric Yang-Mills theory. 

The first calculation of a heavy quark potential in vacuum for N = 4 S YM theory was carried 
out by Maldacena in yfl. There, the potential is obtained from the expectation value of a static 
temporal Wilson loop which, on the gravity side, corresponds to the action of the worldsheet 
spanned by an open string connecting the quark and the antiquark. Due to conformal invariance 
of N = 4 SYM theory, the potential is of Coulomb type: 
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with r is the distance between the quark and the anti-quark, A the 't Hooft coupling and Co = 
r 2 (i)/(2 7r) 3/2 . Soon after H calculations of the heavy quark potential for N = 4 SYM theory 
at finite temperature appeared in J3ll3l. The potential obtained in these calculations starts out at 
small r being close to the vacuum potential of Eq. {T), but rises steeper than the vacuum potential, 
becoming zero at a separation r* = 0.154/nT. For larger separations, i.e., for r > r* ', the authors 
of JilBl argue that the string "melts", and the dominant configuration corresponds to two straight 
strings stretching from the quark and the anti-quark down to the black hole horizon. The resulting 
potential is thus zero for r > r* and has a kink (a discontinuity in its derivative) at r — r*. Here 
I summarize the results presented in fl, where some modifications to the calculations in IHQl 
were proposed in order to avoid the unphysical features of the resulting finite-T potential, namely 
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its trivial infrared behavior plus the presence of discontinuities in its derivative (which leads to 
infinite forces). 

In S U (N c ) gauge theory in Euclidean time, the heavy quark potential or, more precisely, 
the free energy, is calculated through the connected correlator of two Polyakov loops at spatial 
separation r: 
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where j3 = l/T. Eq. (fJJ is the definition of singlet V\(r) and adjoint V c ,dj(r) potentials. In the 
AdS/CFT set up, the color decomposition of Eq. (|2]i can be read off from the different open 
string configurations attached to the quark and the antiquark. One can have a "hanging" string 
connecting the quark and the antiquark or, alternatively, two "straight" strings going from the 
(anti)quark to the black hole horizon, Accounting in AdS indicates that the hanging string con- 
figuration gives V\ (r), since there is just one way to connect the quark and the antiquark, while the 
two straight strings give V a dj(f), since there are N% — 1 different ways in which two independent 
strings can be connected to the N c D3-branes living at the bottom of the AdS space. Moreover, 
the two straight strings can only be connected through graviton exchange in the bulk. Therefore 
the contribution of the two straight string configuration to the potential is 1 fN% suppressed, in 
agreement with the analysis of ]6|] 

We calculate the singlet potential V\{r) according to the definition proposed in [7] in the 
real-time formalism: V\(r) is given by the expectation value of a static (temporal) Wilson loop 
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with the temporal extent of the Wilson loop T — > oo. To calculate the Wilson loop in the 
AdS/CFT set up, we start with the AdSs black hole metric in Minkowski space 
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where die 1 = (dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 , z is the coordinate describing the 5th dimension and L 
is the curvature of the AdS$ space. The horizon of the black hole is located at z = Zh with 
Zh = 1/nT. We want to extremize the open string worldsheet for a string attached to a static 
quark at x l = r/2,x 2 = x 3 — and an anti-quark at x 1 = -r/2,x 2 — x 3 = 0. Parameterizing the 
static string coordinates by X M ■■ 
Goto (NG) action as 



: [X° = t, X x = x, X 2 = 0, X 3 = 0, X 4 = z(x)] we write the Nambu- 
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where z' = dz(x)/dx. Finally, the action in Eq. (0 contains a UV divergence associated to the 
infinite mass of the quarks which has to be subtracted out. Usually, the subtraction contains a 
finite piece as well lHH^l which may be temperature-dependent in the case at hand. Here we 
will use the following subtraction: We define the quark-anti-quark potential by 



V(r) = 



- {S NC (r, T) - Re[S NC (r = oo, T)]} IT. 



(6) 



Our prescription is different from the one used in where the contribution of two straight 

strings were subtracted. It insures that the real part of the potential V(r) goes to zero at infinite 
separations. The definition in Eq. © is consistent with the original prescription outlined in [2] 
to find the heavy quark potential at zero temperature. Importantly, it is also consistent with our 
choice of color representation, since only singlet configurations appear in the definiiton of V(r). 

What follows next is just a problem of classical mechanics: One has to find the string trajec- 
tories that minimize the NG action Eq. (01. In order to do so one has to derive and solve the Euler- 
Lagrange equations with the appropriate boundary conditions, z(x = +r/2) = 0. The details of 
this calculation are explained in [5]. Importantly, the classical solutions can be parametrized in 
terms of the quark- antiquark separation r, the temperature T (both "external" parameters), and 
the maximum of the string along the 5th dimension, z max . The later is determined by the condi- 
tion z' max (r, x = 0) = 0. We obtain the following expression for the heavy quark potential of the 
N = 4 SYM theory at finite temperature: 



V(r) = 
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where F is the hypergeometric function. Importantly, the solutions for z ma x become complex 
for r > r c — 0.870 Zh, leading to complex-valued string trajectories z(x) and potential V(r). 
This led the authors of [j, Ql to abandon their solution for r > r c . We suggest however to 
interpret the complex-valued saddle points as corresponding to quasi-classical configurations in 
the classically forbidden region of string coordinates, as previously suggested in [l]. This is 
similar to the method of complex trajectories used in quasi-classical approximations to quantum 
mechanics. 

Moreover, the solution for Zmax has several branches. In order to select the physical branch we 
impose two conditions: i) The right branch should map onto Maldacena's solution at small r. ii) 
In order to get a sensible the quantum-mechanical time-evolution we require that Im[V(r)] < 0, 
sot that the probability of a state, ~ e Im[£ l ' does not exceed unity. These criteria suffice to single 
just one branch for Zmax- Inserting it into Eq. (0), we readily obtain the real and imaginary parts 
of the resulting potential, plotted in Fig.[T]for two non-zero temperatures, along with the zero-T 
curve for comparison. 




Figure l : The real and imaginary parts of the heavy quark potential plotted as functions of the separation r for several 
different temperatures. We put A = 10. 



Several comments are in order: First, the real part of the potential is a smooth, negative 
definite function, going to zero at large inter-quark separations. This latest feature is a direct 



consequence of our renormalization prescription Eq. If we had used the renormalization 
proposed in HQ], then Re[V(r)] would have gone to a positive constant as r — > oo. Moreover, 
the non-zero temperature curves exhibit a strong screening compared to the zero temperature 
case, while at small r we recover the zero temperature potential of On the other hand, the 
heavy quark potential develops an imaginary part for r > r c . This means the potential becomes 
absorptive, as the qq singlet state may decay in the medium, with an absorption rate that increases 
with the quark-antiquark separation. The existence of an imaginary part in the heavy quark 
potential has been previously observed in perturbation theory in I7J]. 

The nature of the screening exhibited by the lines in Fig.[T]can be better understood from the 
following asymptotic expansions for the real and imaginary part of the potential: 
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Instead of the exponential falloff with r characteristic of Debye screening expected from pertur- 
bation theory and which has been postulated for N — 4 SYM theory at strong coupling in J3l, 
the real part of the heavy quark potential falls off as a power, Re[V(r)] ~ 1/T 3 r 4 , at large r. If 
our hypothesis of using the complex string configurations is confirmed, this would be an inter- 
esting new type of screening for the potential. However, the large negative imaginary part of the 
potential leads to exponential decay with time of the heavy quark pair, which may blur the po- 
tential observable consequences of such new screening mechanism. Also, combining the large- 
and small-r asymptotics we can interpolate the real part of the potential to write an approximate 
formula 
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where the parameter ro given can be interpreted as the screening length. 

Several concluding remarks: First, Debye screening is the result of a perturbative calculation. 
Therefore, it may not be valid for distances of the order of r ~ 1 /(g 2 T), where non-perturbative 
effects may be important. On the other hand, one should also bear in mind that the particle con- 
tent of N — 4 SYM is quite different to that of QCD. In particular it does not include particles in 
the fundamental representation (i.e. quarks), but only in the adjoint (The heavy quarks in our cal- 
culation ought to be regarded as external test charges). Since a charge in the adjoint cannot fully 
screen a fundamental charge, it is reasonable to expect a milder screening. Finally, the power-law 
fall off of the heavy quark potential has also been suggested in perturbative calculations jjj. 
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